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Abstract 



We compute the Chern-Simons current induced by Wilson fermions on a d = 2n + 1 
dimensional lattice, making use of a topological interpretation of the momentum space 
fermion propagator as a map from the torus to the sphere, T d — > S . These mappings 
are shown to fall in different homotopy classes depending on the value of m/r, where m 
is the fermion mass and r is the Wilson coupling. As a result, the induced Chern-Simons 
term changes discontinuously at d+ 1 different values for m, unlike in the continuum. This 
behavior is exactly what is required by the peculiar spectrum found for a recently proposed 
model of chiral lattice fermions as zeromodes bound to a domain wall. 
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Recently a method for simulating chiral fermions on the lattice was proposed by one 
of us [|J. The idea is to implement an odd dimensional d = 2n + 1 theory of Wilson 
fermions with a mass coupling of the fermions to a domain wall. The low energy effective 
theory consists of massless chiral fermions bound to the d — 1 dimensional domain wall, 
without there being "doubler" modes from the Brillouin zone boundary. The anomalous 
Ward identities for these chiral fermions in the presence of background gauge fields can be 
directly measured on a finite lattice; this has already been done numerically for d — 3 and 
the results are in agreement with the continuum anomaly in 1 + 1 dimensions ||. In this 
system, the anomalous divergence of the d— 1 dimensional zeromode flavor currents is due 
to charge flow in the direction normal to the domain wall, even though there is a mass gap 
off the wall — an effect discussed by Callan and Harvey ||] for continuum fermions coupled 
to a domain wall. They pointed out that the Chern-Simons action induced by integrating 
out the heavy fermion modes 0, being proportional to m/\m\, has opposite signs on the 
two sides of the domain wall. This gives rise to a Chern-Simons current in the presence 
of background gauge fields with a nonzero divergence at the domain wall. Furthermore, 
the divergence exactly reproduces the even (d — 1) dimensional anomaly for the single 
chiral fermion zeromode that is bound to the domain wall. In this Letter, we perform the 
Callan-Harvey (CH) analysis for the lattice theory in Euclidian space, where the zeromode 
spectrum is more complicated than in the continuum. 

It is far from obvious that the lattice theory should follow the CH continuum analysis; 
after all, the coefficient of the Chern-Simons action gets 0{1) contributions from arbitrarily 
heavy fermion modes, and the heavy spectrum on the lattice looks nothing like in the 
continuum. In fact, we know the induced Chern-Simons operator must have a coefficient 
very different from the continuum result. While ref. jlj analyzed the spectrum of the 
theory for a Wilson coupling r = 1 and a domain wall height < mo < 2 and found a 
single chiral mode, a recent paper by Jansen and Schmaltz || analyzes the same model 
for general parameters and shows that the spectrum bound to the domain wall changes 
discontinuously with varying mo/r 0. They find that for 2k < \mo/r\ < 2k + 2 — where k 
is an integer in the range < k < (d — 1) — there are (^j^^^j chiral modes bound to the 

domain wall with chirality (— l) k x sign(mo); there are no chiral fermions for |mo/r| > 2d. 
This is quite different than the continuum theory, for which there is a single chiral mode 
for any mo 7^ 0. If the induced Chern-Simons action on the lattice is to correctly account 
for the anomalous divergences of the chiral fermion currents on the domain wall, then 
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evidently its coefficient must also depend discontinuously on mo/r in a very particular 
way. We show in this Letter that that does indeed happen ||. 

The Abelian Chern-Simons action in d = (2n + 1) continuous Euclidian dimensions is 
given by 



— e I H 2n + 1 -r A f) A ■ ■ ■ f) 1 

L CS — e "l-"2n + l / sia 1 Ua 2 ^3 u «2n a 2n + l ' 



(1) 



When a massive fermion is integrated out of the theory it generates a contribution to 
the effective action of the form S e ff = c n Tcs- The coefficient c n can be computed by 
calculating the relevant portion of the graph in fig. [I]. This is true on the lattice as well in 
the weak field, long wavelength limit for the gauge fields. Denoting the fermion propagator 
and photon vertex as S(p) and ik^{p,p') respectively, the graph of fig. [3] yields a value for 
c n which may be expressed as 
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The factor of (n + 1) in equ. is due to the symmetry factor of the graph fig. |T]; the 
factor of i is the product of i n+1 from the photon vertices and (— i) n from relating the 
derivatives in equ. (P to powers of momenta. The p-integration is over the Brillouin zone 
of a hypercubic lattice with lattice spacing a = 1. 

The integral (0) looks very difficult to compute on the lattice for arbitrary n, as both 
S(p) and A p (p, p') are in general rather complicated functions. It is made quite tractable, 
however, by exposing its topological properties. Gauge invariance implies that the photon 
coupling satisfies the Ward identity |§ 



allowing c n in equ.(Ef) to be reexpressed as 
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where the differentiation di is with respect to pj. The fermion propagator may be written 
in the generic form 

S~ 1 (p) = a(p) + ibijp) ■ 7 

= N(p) (cos \9(p)\ +i§(p) -7sin|%)|) (5) 
= N(p)V(p) 

where N(p) = V a 2 + b ■ b, 9(p) = Sarctan(|6|/o), and V(p) is seen to be a 2" x 2 n 
unitary matrix. Provided that S~ 1 (p) doesn't vanish for any p, equ. (ffl) is independent of 
N(p), allowing S~ 1 (p) and S(p) to be replaced everywhere by V(p) and (p) respectively. 
This matrix V(p) is seen to describe a mapping from momentum space — which on the 
hypercubic lattice is the torus T d — onto the sphere S defined by the vector 6{p). The 
homotopy classes of such maps are identified by integers, and so the integral in equ. @) 
has a simple topological interpretation: it is, up to a normalization constant, the winding 
number of the map described by the fermion propagator. 

We now proceed to compute this winding number for the Wilson fermion propagator. 
Although we are ultimately interested in the effective action for lattice fermions in the 
presence of a domain wall, we can compute the effective action far from the mass defect on 
either side by treating the fermion mass as a constant m. Thus we can use the standard 
Wilson propagator 

d d 
S ~ 1 (p) = X/ 7 ^ sin ^ +m + r 1 ^2 [ cos 2V ~ !] • ( 6 ) 

fJ,=l M=l 

Continuous changes in the mass and Wilson coupling, m and r, cannot change the value 
of the winding number except at points for which has a zero for some momentum p. 
Such singular points of the mapping occur only at momenta corresponding to the corners of 
the Brillouin zone, and then only for m/r = 0, 2, • • • , 2d. The Chern-Simons coefficient c n 
as a function of m/r must therefore be piecewise constant, changing only at these critical 
values. Furthermore, for fixed r, V(p) — > ±1 as m — ► ±oo, so we may deduce that 

c n (m, r) = for m/r < , m/r > 2d . (7) 

To compute c n for < m/r < 2d, we need only evaluate the derivative of the integral 
in equ. with respect to m across the critical values m/r = 0, 2 ... , 2d. This task is 
simplified by the fact that c n (m + dm, r) is unchanged as one deforms dm in a p-dependent 
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way so that dm(p) vanishes for all p except in the vicinity of the Brillouin zone corners; 
these points are denoted by p = £a k \ the a = 1, . . . , (^j vectors with k nonvanishing 
components equal to n. Therefore we need only evaluate the integrals in infinitesimal 
regions near the Brillouin zone corners £ a {k) . After some algebra this yields 



dc n (m,r) _ i(-l) n 2 n fc _d f d 2n+1 p (m - 2rk) 

~ (n + 1) ff~ d^ J dn (2tt) 2 -+ 1 [\p\2 + ( m _ 2rk) 2 ] n+1 



(8) 
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where the integration dfi is over the infinitesimal region \p\ < e — > 0. This may be trivially 
integrated with respect to m, given the boundary values (^), to yield the Chern-Simons 
coefficient for Wilson fermions: 

c n (m,r) = —, — - }) : — ry^(-l) fc if I 7— — — tt • (Wilson fermions) (9) 
v ' ; 2(n + l)(2n) n nl ^ ' \kj\m-2rk\ K 1 y> 

This is to be compared to the continuum result for a fermion of mass m, computed by 
inserting the continuum free propagator in equ. (|j) and integrating over continuum mo- 
mentum space: 

i( — l) n TTl 

c n (m) = — — — . (continuum result) (10) 

2{n + l)(27r) n n! \m\ 

Our computation is exact in the limit of small and adiabatic external gauge fields, even 
for m of 0(1) in lattice units; it is readily generalized to non-Abelian gauge fields. 

So far we have only considered lattice fermions with a constant mass; our primary 
interest though is in studying current flow in the presence of a domain wall mass term for 
the fermion. Following Callan and Harvey, we assume that sufficiently far away from the 
domain wall on either side, we can treat the mass as constant and compute the induced 
current by varying the Chern-Simons action, 

jCS = SS eff (n + 1) 
A* 2 n ' )t ^ a l a 2- a 2n + l ra l a 2' a 2n-l Q 2n ' ^ ' 

where the factor of 2 n arises from replacing dA by the field strength F. Since c n depends 
on the mass m(s) which is space dependent, J^ CS has a nonzero divergence corresponding 
to current flow normal to the domain wall. In their continuum analysis, Callan and Harvey 
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pointed out that with c n given by equ. (|10|), there was equal current flow toward the wall 
from each side that exactly compensated for the anomalous divergence of the chiral fermion 
current along the wall's surface. For lattice fermions in the presence of the domain wall, we 
must use the value (||) for c n (m, r). Thus on the side of the domain wall for which m(s)/r is 
negative, the Chern-Simons current vanishes. On the other side, where m(s)/r — > |rao/r|, 
the current either vanishes if \mo/r\ > 2(d + 1), or else, for 2£ < \mo/r\ < 2(£+ 1), the 
current is given by 



.-(lattice) £ ( _ 1)t (*)_5W<) 



JCS (continuum) 

v v ' k=0 k=t+l 



(12) 



2(-l 



d-1 



For example, for I = (a domain wall height satisfying < \mo/r\ < 2) the lattice current 
vanishes on one side of the wall, while having twice the continuum magnitude on the other 
side; in this case the total divergence of the continuum and lattice currents are the same. In 
general the lattice Chern-Simons current vanishes on the side of the domain wall for which 
m(s)/r is negative, while on the other side has the correct magnitude to compensate 
for the anomaly of not one positive chirality zeromode at the domain wall, but rather 
(^£^J zeromodes with chirality (— l) e x sign(mo/r). This result agrees precisely with 
the zeromode spectrum bound to the lattice domain wall found by Jansen and Schmaltz 
0. It also agrees with a numerical computation of the Chern-Simons current on a d = 3 
finite lattice that we have performed (see ref. 0), which exhibits the behaviour peculiar 
to the lattice that the Chern-Simons current flows on only one side of the domain wall, as 
well as the discontinuous dependence of its magnitude on mo/r. 

Our analysis confirms that the model of ref. [l[ correctly reproduces the continuum 
anomaly for chiral fermions on a finite lattice in the presence of weak gauge fields, even 
for a domain wall height mo = 0(1) in lattice units. 
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Figure Captions 

Fig. 1. The Feynman diagram in 2n + 1 dimensions contributing to the induced Chern- 
Simons action for Abelian gauge fields; Y^i=\ q% = 0. Graphs with multiple 
photon vertices peculiar to the lattice do not contribute, having the wrong Lorentz 
structure. 
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